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,— i Abstract. We obtain the mean field limit and the propagation of chaos for a system of 

particles interacting with a singular interaction force of the type l/|x| a , with a < 1 in 

dimension d > 3. We also provides results for forces with singularity up to a < d — 1 but 

qj with large enough cut-off. This last result thus almost includes the most interesting case of 

Ph Coulombian or gravitationnal interaction. 
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1 Introduction 

The N particles system. The starting point is the classical Newton dynamics for point- 



particles. We denote by (Xi, . . . , Xn) the position of the particles in R . , and by (Vi, . . . , Vj 



A' J 



their velocities in IR d . Assuming that particles interact two by two with the interaction kernel 

Fix), one finds the usual 
> 

8 v i =E N {x i ) = Y,^nx i -x j ). (i.i) Egjjg 

O 

We use the so-called mean-field scaling which consist in keeping the total mass (or charge) of 

order 1 thus formally enabling us to pass to the limit. This explains the 1/N factor in front 

.1^ of the force terms. This implies corresponding rescaling in pom tipn, velocity and time. 

There are many examples of physical systems following (EnTJT^The best known concerns 



Coulombian force F(x) = C x/\x\ d ~ l , which serves as a guiding example and reference. Those 
describe a plasma, or for C < gravitational interactions, in which case the system under 
study may be a galaxy, a cloud of star or galaxies (and thus particles can be "stars" or even 
"galaxies"). For simplicity we consider here only a basic form for the interaction. However 
the same techniques would apply to more complex models, for instance with several species 
(electons and ions in a plasma), 3-particles (or more) interactions, models where, the force 
depends also on the speed has in swarming models like the Cucker-Smale one [CDPJOSaj... 
For convenience, we also use the notation Z, = (Aj, Vi) for the solution and 
Z° = (X°, V?, . . . , A°, 1/°) for the given initial conditions. 
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Finally let us mention that sometimes the kernel F in fact depend on the number of particles. 
This might seem quite strange from the physical point of view but is in fact very common 
for numerical simulations (to regularize the interactions). When such a cut-off is used, we 
will denote by F^ (some N dependent mollification of F) the force field. 

Or 



The Jeans-Vlasov equation. At first glance, the system (llTtrinight seem quite reason- 



able However many problems arise when one tries to use it for practical applications. But 
in our case, the main issue is the number of particles, i.e. the dimension of the system. For 
example a plasma or a galaxy usually contains a very large number of "particles", typically 
from 10 10 to 10 23 , which makes solving (TTTrunpossible. . n nF 

As usual in this kind of situation, one would like to replace the discrete system (ITT pay a 
continuous model. In that case this model is posed in the phase space, i.e. it involves the 
distribution function f(t,x,v) in time, position and velocity. The evolution of that function 
f(t, x, v) is given by the Jeans-Vlasov equation 



d t f + vVJ + E(x)-V v f = 0, 
E ( x ) = f*d p(t, y) F(x - y) dy, 
P(t' x ) = f v f(t,x,v)dv, 

where here p is the spatial density and the initial density f° is given. 

Our whole concern in this art icle j &J>o understand when and in which sense, 

be seen as a limit of system (ITTr^rhis question is of obvious importance for theoretical 



1.2) eq:vlasov 



Eq. dO 



vlasov 



) can 



reasons, to justify the validity of the Vlasov equation for example. 

It also plays a role for numerical simulation, and especially Particles in Cells methods which 
introduce a large number (roughly around 10 6 or 10 8 , to compare with the order 10 10 to 
10 23 mentioned above) of "virtual" particles in order to obtain a particles system solvable 
numerically. The problem in that case is to explain why it is possible to correctly approximate 
the system bou sing much fewer particles. This would of course be ensured by the convergence 
of (O) to Of. 



ea-^ylasov 



Formal derivation of Eq,(lT]2[ from (ITTTrT One of the simplest way to understand 
formally how to derive Eq. (ITzj is to introduce the empirical measure 



1 N 



(t),Vi(t)- 



Then if J X^ 1 Vj) is solution to (ITTFrand if there is no self-interaction : F(0) = 0, then p z N 



solves (ITTz) in the sense of distribution. Formally one may then expect that any limit of p z N 
still satisfies the same equation. 



The question of convergence. The ^previous formal argument suggests a first way of 
rigorously deriving the Vlasov equation (fO). Take a sequence of initial conditions Z° (to be 
given for every number A^ or a sequence of such numbers) and assume that the corresponding 
empirical measures converges (in the usual weak-* topology) 



^(0)— >f(x,v) 



One would then try to pro veth at the empirical measures at later times /J*%(t) weakly converges 
to a solution f(t,x,v) to (pTTtZ) with initial data f°. 



In other words, is the following diagram commutative? 

/4(o)^/(o) 



Npart 



VP 



Note that for singular kernels F (among which Coulombian and gravitationnal interactions), 
one does not expect to be able to do that for any initial conditions. First of all solutions to 

je a _; ,vl_asov _„ . _ -„ 



(1172) do not exist in general if the initial data / is only a measure. And even if / is smooth 



but a small amount of the particles are initially concentrated in a small region, problems will 
likely occur (as the interactions blow up) 

Of course it is not obvious what are admissible initial conditions or how to precise that. The 
simplest way is to give some properties that the initial conditions must satisfy. For instance 
one could ask the initial positions and velocities to be uniformly distributed (on a grid for 
numerical simulations). However from the point of view of statistical physics, this is too 
restricitve. This leads us to the notion of propagation of chaos. 

Propagation of chaos. In most physical settings, one expects the initial positions and 
velocities to be selected randomly and typically independently. In that case the law is initially 
given by (y )®^ (i.e. randomly and independently with profile /°). Note that the empirical 
measure at time is then close to /° with large probability (in some weak norm, See the 
Proposition pT for a precise version). 

The propagation of chaos was formalized by Kac's in l [Kaco6J and goes back to Boltzmann 
and its "Stosszahl ansatz" . A standard reference is the famous famous course by Sznitman 

Szrn t.nia n 

fSznylJT It is mainly use in probabilistic system, where some randomness is introduced in 
the dynamics of the particle, but is also relevant in our context. m n nF 

Denoting by f N (t,x±,Vi, . . . ,xn,vn) the image by the dynamics (ITT Pot the initial law 
(f°)® N , one may define the fc-marginals 

cNu _ _ \ / fN 



L f 



f k (t, Zi, . . . , z k ) = f (t,zi,...,z N )dz k+1 ...dz N . 

J R 2d(N-k) 

According to the general definition, the propagation of chaos means that for any fixed k, 

fk(t) wm converge weakly to (f)® k as iV — > oo. In fact it is sufficient that the convergence 

holds for only one k > 2. 

As Sznitman shows, it is also equivalent to say that the empirical measu ^j^ii(|) converge in 

law towards the constant variable /(£). In fact we will give in theorem Iff and E a quantified 

version of the convergence in probability of (i%(t) towards /(£). 

Shortly, this is possible because, the marginals can be recovered from the expectations of 

moments of the empirical measure 

f£ = E(/4(£,Xi,Wi) . . .fj, z N (t,x k ,v K )) + O ( — 

a result sometimes called Grunbaum lemma. For detailed explanations about quantification 
of the equivalence between convergence of the f k and the convergence of the law of the 



empirical distributions, we refer to |HMi2j7 This quantified equivalence was for instance 

used in the recent and important work of Mischler and Mouhot about Kac's program in 

kinetic theory [MM 11 J: 

Of course we do not expect prop agat ion of chaos to hold for any initial distribution /°, for 

instance if f° 

f° G L°°. 

for System (111 .11) with a singular force kernel, namely whether one can even expect to have 



is too singular for (1172) to have a solution. Hence we limit ourselves here to 



Crra 



Well posedness Tor System (OT. We have not mentioned yet the most basic question 



solutions to the system for a fixed number of particles. Indeed, because of the singularity, 

the usual Cauchy-Lipschitz theory cannot be applied. 

First of all for the type of singularity that we will handle here, the answer is relatively easy 

as it would be simple to show that velocities remain bounded and that particles cannot 

collide for almost all initial data. We do not give a specific proof as that result is a simple 

consequence of our analysis. 

However in more singular cases (especially repulsive ones), the question remains. The classical 

approach to this problem uses the so-called theory, of renormalized solutions developed by 

DiPerna-Lions. We refer to IJDL89Jand to Hauray |Hau04] for this specific problem. 



Previous results with cut-off or for smooth interactions. The convergence and the 
propagation of chaos are known to hold for smooth interaction forces (F G C 1 in ,£gner 
or at least W 1 ' 00 ) since ^he end of the sevent^es^nd the works of Braun and Hepp flBHVV 
Neunzert and Wick NW80J and Dobrushin |Dob79j. Those articles introduces the ] 
and the formalism behind mean field limits, we also refer to the nice book by Spohn 
Their proofs however rely on Gronwall iyp axstimates and are connected to the fact that 
Gronwall estimates are actually true for (IHTruniformly in N if F G W 1 ' 00 . Those makes them 




impossible to generalize to any case where F is singular (including Coulombian interactions 
and many other physically interesting models). 

Instead, by keeping the same general approach, it is possible to deal with singular interactions 
with cut-off. For instance for Coulombian interactions, one could consider 



F N (x) = C 



x 



(\x\+e(N)y 



or other type t oi s ecularization at the scale e(N). 

The system (fnTrcToes not have much physical meaning but the corresponding studies are 



crucial to understand the convergence of numerical methods. 

lyTrvtory 

For particles initially on a rnesh. we refer to the works of Ganguly and Victory |GV89[, 



Wollman flWolOO' 
cut-off) 



and Batt 



BatOlj (the later gives a simpler proof, but valid only for larger 
Unfortunately they had to impose that limAr_ > . 00 £(./V)/iV~ 1 / o! = +00, meaning that 



the cut-off for convergence results is usually larger than the one used in practical numerical 
simulations . Note that the scale N~ l / d is the average distance between two neighboring 
particles in position. 

Of course propagation of chaos cannot be proved in those cases as the particles are on a 
mesh initially and hence cannot be taken randomly. Moreover, we emphasize that the two 
problems with particles initially on a mesh, or with particles not equally distributed seems 
to be very different. In the last case, the previously mentioned results do not apply, and 



IVirrtorvp 

Ganguly, Lee and Victory |GLV9'1J are only able to prove the convergence for a very large 
cut-off e^ (In N)-\ 



Previous results for 2d Euler or other macroscopic equations. A well known case, 
very similar at first sight with the questiomhere, is the vortices system for the 2d incom- 
pressible Euler equation. One replaces (111. if by 



Xi = j^^2®i a, F(Xi - Xj) 

3+i 



where F is still the Coulombian kernel (in 2 dimensions here) and ai 
system to converge to the Euler equation in vorticity formulation 

d t uj + div (u uj) = 0, div u — 0, curl u = 



;i.3) 



±1. One expects this 



U. 



The same questions of convergence and propagation of chaos can be asked in this setting. 
Two results without regularizatipn for the true kernel are already known. The work of 
Goodman, Hou and Lowengrub, |GHL90j and |GH91[, has a numericaLpnint of view but use 
the true singular kernel in a interesting way. The work of Schochet fSch96] uses the weak 
formulation of Delort of the Euler equation and provethat empirical measure with bounded 
energy converges towards weak measures solution to (DPtrr Unfortunately, the possible lack of 



uniqueness of Euler equation in the class of measures do not allows to deduce the propagation 
of chaos. _ , , 

As equations like (OEFare notoriously harder to deal with than kinetic equations like (IT72F7 



vortex 



(1.4) | Euler 



vlasov 



one could expect similar results for our problem. Unfortunately, the mean field limits are 
more difficultin the phase space. Thereare several reasons for that, in particular the fact 
that system ( |ll.'l[ ) is second order while (EnSTi^ first order. This implies Jdiat collisions or near 
collisions (in physical space) between particles are very common for ( li.'lp even for repulsive 
interactions and rare for (lortat least for vortices of the same sign). 



For example, the references mentioned above use the symmetry of the forces in the vortex 
case, a symmetry which does not exist in our kinetic problem. The force is still symmetric 
with respect to the space variable, but there is now a velocity variable (second order again) 
which breaks the argument used in the vortices case. For a more complete description of the 
vortices system, we refer to the references already quoted or to |Hau09j , which introduces in 
that case techniques somewhat similar to the one used here. 



Previous, results in singular cases without cut-off. Let us first mention that the 

>?q :,vIasov _ 

ngular, 
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equation (1172) at the limit is now well understood, even when the interaction F 



iBol 



eluding the Coulombian case. The existence of weak solutions .goe^ back to |DabVSll or 
ArsVoj . Existence of global classical solutions is proved in Pfa92j, ISchylJ (see also |Hor93j) 
and at the same time in I JLP91J . Of course those results require some smoothness on the 
initial data f° (for instance compactly support and bounded). 

To our knowledge however, the only mean field limit result available up to now is fHJOVJT 
This proves the convergence (not the propagation of chaos) provided that 



07 



The interaction kernel F behaves like 



x\ 



~ a with a < 1. 



• The particles are initially well distributed, meaning that their minimal distance in 
phase space is of the same order as the average distance between neighboring particles 

The second assumption is all right for numerical purposes but does not allow to consider 
physically realistic initial conditions (as per the propagation of chaos property). This as- 
sumption is indeed not generic for empirical measures chosen with law (f°^® N [i.e.. it is 
satisfied with probability going to in the large N limit). 



thm:prob 



Our result without cut-off. In the present article, we keep the same conditions on the 
inteaction kernel, but require only a much weaker assumption on the minimal distance be- 
tween particles. This allows us to prove the propagation of chaos, for forces satisfying a 
(S' a )-condition : 

(S a ) V ie i rf , \ F (x)\<£-, |VF(x)| < -£tt, (1.5) 



\x\ 



\x\ 



with a < 1. Our precise result without cut-off is the following 

Theorem 1. Assume that d > 3 and that F satisfies a (S a )- condition with a < 1 



Choose 
vu initial condition / u G L°° with compact support and totaLnmss one for the Vlasov equation 
(ETT2). For each N E W , consider at the particles system ( |ll.l[ ) with initial positions (X«, Vi)i<N 
chosen randomly according to the probability (J )®^. Then for all T > 0, all 



2 + 2a 



< 7 < 1 and < s < 



7<J - (2 - 7)0 - 2 



a) 



d + a ' 2(1 

there exists positive constants Co(f, F), C s (7, s, f, F) such that for N large enough (InN >T) 



F(3te[0,T], W 1 ( m (t)J(t))> 



3e 



C t 



JV7/(2d) J 



< 



9jl 

N s 



FH 



vlasov 



where f(t) is the unique strong solution of the Vlasov equation (fffz) with initial condition 
f°, (the constant C s blows up when s approaches its maximum value) and W\ denotes the 1 
Monge-Kantorovitch- Wasserstein distance. 

The notation C(f) means that the constant depends on the function / (essentially via con- 
served quantities like ||/||oo and also the size of its support) on the whole interval of time 
under consideration, here [0, T]. 
The conditions on 7 and s are not completely obvious, but it can be checked that if a < 1 



< 1 so that admissible 7 exists. And for an admissible 7, the quantity 



2+2a 
d+a 

is also positive, so that admissible s also exists. , 

* im:prob 



and d > 3, 

7<f-(2-7)q-2 

2 ( 1 + Q ) -■ ,---- 

Roughly speaking, under the assumption of Theorem [TJ the probability of finding a deviation 
strictly larger than the average inter-particle distance N~ l l 2d is small. 

Remark 1. Unfortunately, we are not able to provide a similar 



@. . result Jor d = 2, even if a 
im:proP ' J 

is empty in that case. 



Even if the theorem is stated probabilistically in terms of 1 
of the proof is a deterministic theorem (See the Theorem B 
which has generic assumptions with respect to the law (f°)® J 



propagation of chaos 
uiKTlerer 



the core 
stated in the second section) 



Thanks to the deterministic 
result we can also construct explicit sequences of initial conditions for which the convergence 
towards the Vlasov equation will holds (for instance, particles well choosen on a mesh, but 
not only). 



eq:Calpha 



1.6) eq:thmprob 
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The improvements with respect to |HJ07J7 The major improvement is of course the 
much weaker condition on the inital distribution of positions and velocities. We are hence 
able to show the propagation of chaos, which is again the crucial property for applications 
to physics. „ . , ._ 

We managed to simplify the proof with respect to the previous article |HJ0V|7~considerably 
so in the long time case which was quite intricate before and does not require any special 
treatment here. 

Finally our new result is almost quantitative. For large enough N, it actually tells how close, 
in Wasserstein distance, the empirical measure is from the limit. In fact if one does not use 
random initial conditions (and hence we do not need the large deviation result mentioned 
above), Theorem p\ gives a precise rate of convergence, which is of course quite useful from 
the point of view of numerical analysis. 

Unfortunately, the condition on the kernel F is still the same and does not allow to treat 
Coulombian interactions. There are some physical reasons for this condition: for instance if 
F = — V0 then it guarantees that the potential is bounded. We refer to |BHJ10J for some 
ideas in how to go beyond this threshold. 



probcutof f 



The result with cut-off. The result„,with cut-off presented here is in one sense slighty 
weaker than the previously known result |GV89f, since we just miss the critical case a = d — 1. 
But it has also some major advantages, especially if we are not interested by numerical 
simulation. First of all, it is valid for random initial configurations with propagation of chaos 
and not only for well distributed initial positions and velocities (on a mesh). Secondly, for 
a larger but close to one it is valid for smaller cut-offs, much smaller than average (and also 
minimal) distance between particles. 
The result is stated for forces depending on N and satisfying the following condition 



Wm. 



It) 

Hi] 



F satisfy a (S a ) — condition 
V|x| >N~ m , F N (x) = F(x) 
V|x| < N~ m , \F N {x)\ < N~ ma , 



which essentially means that the interaction kernel is regularized at scales lower than N~ m . 
Note that in fact we would not need any estimate on the gradient of Fn for very small x. 
The result would still be true if Fn only converges to F for large enough x, with an error 
satisfying ||Fjv — F||i < N~ 1 / 2d . The proof could be adapted to that case, but for simplicity 
we choose this presentation. We also point out that one would like to take m as large as 
possible if we want to be close to the dynamics without cut-off. 

Theorem 2. Assume that d > 3, 7 G (0, 1) and that Fn satisfies a (S^)- condition for some 
1 < a < d — 1 with a cut-off order m such that 



m < m 



1_ 
2d 



mm 



d - 2 2d - 1 



a 



o 



Choose anymitial condition f° G L°° with compact support and totcdmass one for the Vlasov 



equation (1X72). For each N G N* ; look at the particles system (\T71 Jwith initial positions 
(Xi, Vi) n <N chosen randomly according to the probability (J )®^. Then for all T there exists 
positive constants Co(f, F), 61(7, m, f, F) and 02(f) such that for N large enough (InN >T) 



F 



4gCo* \ 



3tG[0,T], Wl Mt),f(t))> 1 ^- d ) 



< CoN" 1 e 



-C 1 N > 



1.7) eq:Ckappa 



where, A = min (l — 7, ^jf) and f(t) is the unique strong solution of the Vlasov equation. 



^ 



(lO) with initial condition f° 



In dimension d = 3, the minimal cut-off is given by the order of m* = ^ min ((a — 1) _1 , 5a _1 ). 
As 7 can be chosen very close to one, for a larger but close to one, the previous bound tells us 
that we can choose cut-off of order almost iV _5//6 , i.e. much smaller than the likely minimal 
inter-particles distance in position space ( of order iV _2//3 , see the third section). 
With such a small cut-off, one could hope that it is almost never used when we calculate 
the interaction forces between particles. Only a negligible number of particles will become 
so close to each others during the time T. This suggests that there should be some way to 
extend the result of convergence without cut-off at least to some a > 1. 
Unfortunately, we do not know how to make rigorous the previous probabilistic argument 
on the close encounters. First it is highly difficult to translate for particles system that are 
highly correlated. To state it properly we need infinite bounds on the 2 particles marginal. 
But obtaining such a bound for singular interaction seems difficult. Moreover, it remains 
to neglect the influence of particles that have had a close encounters (its trajectory after a 
encounter is not well controlled) on the other particles. 

Let us also mention that astro-physicists doing gravitational simulations (a = d — 1) with 
tree codes usually use small cut-off parameters, lower than JV _1 / d by some order. See fBehOOj 
for a physical oriented discussion about the optimal length of this parameter. 

|ea-Lvlasov 



A short sketch of the proof. As mentioned above, the Vlasov equation (tLzj is satisfied 



■A 



by the empirical distribution /ijy of the interacting particle system provided that F(0) is set 
to 0. Hence the problem of convergence can be reformulated into a problem of stability of 
the empirical measures - seen initially as a measure valued perturbation of the smooth profile 
f° - around the solution fit) of the Vlasov equation. 
Our proof uses 3 ingredients to obtain this stability 

• Show that with large probability the empirical measure at the initial time /x^-(O) is very 
close to f° and that the particles are not too badly distributed in the phase space. 



o-ylasov _ Jeo-: 



Compare the solution f(t) to (lO) with f° as initial data to the solution /jy to (Of 



vlasov 



with a regularization of fi%(0) as initial data ((Dirac masses are replaced by "blobs"). 

• Control the distance (in some appropriate Monge-Kantorovitch-Wasserstein distance) 
between n z N {t) and /jv(i)- 

The first two steps are not overly complicated because they rely on previous known results 
(the second is a standard stability result for Vlasov-Poisson for example) and rather simple 
probabilistic estimates. The difficulties are hence concentrated in the third step- 
This uses a deterministic result that has to be more precise than the one in |H JOVJT First of 
all it has to allow for a much smaller minimal distance in phase space between the particles 
at the initial time (and thus as well at later times). Notice that in some average sense the 
particles are still reasonably well distributed: For instance, the regularization /at(0) of /i^-(O) 
at a small scale is bounded in L°° . 

The other complication comes from that we want an explicit bound on the distance between 
/jv(£) and n z N {t) in the W^ distance as we will need it to control the distribution of particles 
at time t. This is necessary to control the difference between the force terms and hence the 
evolution of the distance between /jv(^) and ^ z N {t). 



We remark that the use of the infinite MKW distance is important. We were not able to 
perform it with other MKW distance of order p < +00. It may seem strange to propagate a 
stronger norm for a problem with low regularity but in fact it turns out to be the only MKW 
distance with which we can handle a localized singularity in the force and Dirac masses in 
the distribution. 

In essence controlling the distance between /jv(£) and n z N {t) requires to prove that the dif- 
ference between the force terms acting on /jy and on p,%{t) is bounded by this distance plus 
a small correction; then one may conclude by a Gronwall-like argument. There are several 
major complications in the proof. 

First of all one has to deal with large oscillations in the discrete force term. This is due to 
particles that may come very close one to another in physical space or even collide, while 
remaining at a reasonable distance in phase space. The force term may become very large in 
that case but only for a very short time as the particles do not remain close. To solve this 
problem, we have to average those oscillations over a short time interval. 
Then when one compares the force terms it is necessary to distinguish between 3 domains of 
interacting particles: 

• Interaction between particles far enough in the physical space and that remain far 
enough over the short time interval where we average. This is the simplest case as one 
does not see the discrete nature of the problem at that level. The estimates need to 
be adapted to the distance used here but are otherwise very similar in spirit to the 
continuous problem or other previous works for mean field limits. 

• Interaction between particles close enough in the physical space, far enough in the 
phase space over the short time interval where we average. Here we start to see the 
discrete level of the problem and in fact we cannot compare anymore the force term 
for the continuous Vlasov equation with the one for the particles' system. Instead we 
just show that both are small. This is simple for the continuous term as the size of 
the domain is itself quite small. For the discrete force term, it is considerably more 
complicated. First of all one has to bound the average over the short time interval of 
the interaction term between any two particles in that situation. Roughly speaking if 
those have a relative velocity of order v one expects the interaction to behave like 



t+T 



ds 



\5 + sv\ a 

where r is the size of the average in time and S is the minimum distance between the 
two particles over [t, t + r] . This is where the condition a < 1 first comes into play 
as it allows to bound the previous integral independently of S, provided that v is large 
enough. To conclude this part it is finally n^ce^sary to sum the previous bound over all 
particles in the corresponding domain. In |H JOVJ, we could directly bound the number 
of particles there; here as we do not have the same lower bound on the minimal distance 
in phase space, we have to use the distance between f N and \x z N (which is why we need 
the Wnr, distance). 



00 



Interaction between particles close enough in phase space over the short time interval 
where we average. In |H JOVJ this case was relatively simple as the diameter of the 
corresponding domain in phase space was of the same order as the lower bound on 
the minimal distance (still in phase space) that we were propagating. Here this lower 

9 
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bound is much smaller (of the order of the square of the, diameter) . This is where the 
main technical improvement lies with respect to l |HJ0V| . 

Organization of the paper. In the next section, we introduce the notations, and state 
the deterministic results on which the propagation of chaos relies. In the third section, we 
explain how to obtain the propagation of chaos from the deterministic results. The fourth 
section is devoted to the proof of the two deterministic theorems. 

2 Notations and other important theorems 

2.1 Notations and useful results 

We first need to introduce some notations and to define different quantities in order to state 
the result. 

• Empirical distribution /j.^ and minimal inter-particle distance d^ 

Given a configuration Z = (X i? V})j<jv of the particles in the phase space M MAr , the associated 
empirical distribution is the measure 



(*N 



^I> 



Xi,V r 



An important remark is that if (Xi(t,),Vi(t))i<N is a solutiomof the system of ODE (IT7TJ7 
then the measure fJ.%(t) is a solution of the Vlasov equation (ET2), provided that the inter- 
action force satisfies F(0) = 0. This condition is necessary to avoid self-interaction of Dirac 
masses. It means that the interaction force is defined everywhere, but d^continuous and has 



ODE 



a singularity at 0. In that conditions, the previously known results |BH77]7|NW80j cannot 
be applied. 

For every empirical measure, we define the minimal distance d% between particles in phase- 
space: 

d z N {ti N ) = minflXi -Xj\ + \V t - Vj\). (2.1) 

This is a non physical quantity, but it is crucial to control the possible concentrations of 

particles and we will need to bound that quantity from below. 

In the following we often omit the Z superscript Before going on, in order to keep "simple" 

notations. 

• Infinite MKW distance 

First, we use many times the Monge-Kantorovitch-Wasserstein distance of order one and 
infinite. The order one distance, denoted by ||(t,k classical and we refer to the very clear 
book of Villani for definition and properties 
widely used, so we recall its definition : 



eq:dmin 



The second one denoted W^ is not 



Definition 1. For two probability measures \i and u on X, a polish space, with II(/i, v) the 
set of transference plan from [i to v. 

Woo(/i, v) = inf{A — esssup \x — y\ \\ G II}. 
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In one of the few works on the subject |CDPJ08bf Champion, De Pascale and Juutineen prove 
that if \i is absolutely continuous with respect to the Lebesgue measure £, then at least one 
optimal transference plane is given by a optimal transport map. In other words there exists 
a measurable map T : Supp/i — > X such that (Id,T)#fi e IT (it implies in particular that 
T#fi = v) and 

Woo^i, v) — fi — esssupj, \Tx — x\. 

Although that is not mandatory (we could also work with optimal transference planes), we 
will use this result that will greatly simplify the proof. 

Optimal transport is useful to compare the discrete sum of the N particles dynamics to the 
integrals of the continuous Vlasov system. For instance, if / is a continuous distribution and 
/ijv & n empirical distribution we may rewrite the interaction force of hn using a transport 
map T = (T X ,T V ) of / onto /j,n 



iyjF(.Y«-A-») = /F(Al 

i+3 



T x {y,w))f(y,w)dydw. 



Note that in the equality above, the function F is singular at x — 0. Using infinite MKW 
distance, the singularity is still localized "in a ball" after the transport. The term under 
the integral in the right-hand-side has no singularity out of a ball of radius W / oc (/, vn) in x. 
Others MKV distance of order p < +00 destroys that simple localization after the transport, 
which is why it seems more difficult to use them. 

• The scale s. We also introduce a scale 



e(N) = N~^ 2d . 



(2.2) 



for some 7 G (0, 1) to be fixed later but close enough from 1. Remark that this scale is larger 
than the average distance between a particle and its closest neighbor, which is of order N~ l l 2d . 
We shall do a wide use of that scale in the sequel, and will often define quantities directly in 
term of e rather than N. For instance, the cut-off order m used in the (S^)-condition may 
be rewritten in term of s, with m := — m. 



eps 



i) F satisfy a (S a ) — condition 
(S£) ii) V\x\>e™,F n (x) = F(x) 

Hi) V|x| <e m , \F N (x)\ <e-™ a , 



(2.3) eq:Ckappa J 



L>2</ 



• The solution fa of Vlasov equation with blob initial condition. 

Now we defined a smoothing of fi^ at the scale e(N). For this, we choose a kernel <p : K"" 1 — > 
with compact support in [— |, ^] 2d and total mass one, and denote e (-) = e~ 2d (f)(-/e). The 
precise choice of 4> is not very relevant, and the simplest one is maybe <j> = l r _i iud- We use 
this to smooth n^ and define 



I Il2d. 
2'2< 



./. 



N 



Hn * 



>e(N), 



Jeq : . ylasoy 

and denote by fa{t,x,v) the solution to the Vlasov Eq. (IO) for the initial condition f 



f0 

in- 



The interest of f/v ,is that we may assume that it belongs to L°° (see the deviations estimates 

J J pf.op riargedev 7 1 

of the Proposition Hm the Appendix). It allows to use standard stability estimates to control 



its W\ distance to another solution of the Vlasov equation (See Loeper result |Loc06j). 



(2.4) eq:deffN 
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thm: deter 



thm: cutoff 



2.2 Statement of the deterministic result without cut-off 

S uiKprob 
the 

randomness comes only from the choice of the starting initial data. Precisely, the probability 

on the initial conditions is used to ensure that some conditions on minimal inter-particle 

distances and MKV distances are satisfied with large probability But, once that conditions 

are fulfilled, we are able to propagate them with deterministic estimates. 

The following theorem shows that the particles system may be approximated by the solution 

of the Vlasov equation with the "blob" distribution fy as initial conditions, provided that 

two conditions on the minimal inter-particle distance d/v(0) an d the infinite norm of /jy are 

satisfied. 

Theorem 3. Assume that d > 2 and that the interaction force F satisfies a {S a ) condition, 
for some a < 1 and let < 7 < 1. Assume also that the initial empirical distribution fi° N of 
the particles and its e- enlargement f% satisfy : 



1) d° N := d N (fj, N (0)) > e 1+r = JV-T( 1+r )/ M for some r G (l,r*) where r* := 
H) ||/at||oo < Coo, a constant independent of N, 



d-l 

l+a> 



Hi) For some R > 0, ViV G N, Supp/x^ C B(0,R), the ball of radius R and center of 
R 2d . 

Then for any T > 0, there exists two constants Co(R, C^, F, T) and Ci(R, C^, F, 7, r, T) 
such that for N > e ClT the following estimate is true 



V£G[0, T], Woo(M*),M*))< 



pCot 



Remarks. This is a inequality of the type W^t) < W / oo (0)e c ' t , where the value of ^00(0) 
has been bounded by JV~ 7 / 2o! . But that last bound is true since /$• is a blob approximation 
of fi° N , with blob contained in balls of radius N~ 1 ^ 2d around the Dirac of jx q n . 
The previous theorem is valid in dimension 2. But unfortunately, its conditions are not 
generic in that case if the initial conditions are choosen independantly. This is why we 
cannot conclude to propagation of chaos for d = 2. 

2.3 Statement of the deterministic result with cut-off 

■ fchm : probcutof f 

As in the case without cut-off, the probabilistic result |2J relies on a deterministic result, much 
simpler with cut-off since it does not need any control on the minimal inter-particles distance. 
The result is the following 

Theorem 4. Assume that d > 2 and that the interaction force F = Fn satisfies a (S^), for 
some 1 < a < d — 1, with a cut-off order satisfying 

7 . (d-1 2d-l 

m < m := — : mm , 

2d \a — 1 a 

Assume also that the initial empirical distribution of the particles n° N and its e enlargement 
f N satisfy : 



(2.5) eqrthml 
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i) ||/at||oo < Coo, a, constant independent of N, 

ii) For some R > 0, ViV G N, Supp//jy C B(0,R), the ball of radius R and center of 
R 2d . 

Then for any T > 0, there exists two constants Cq{R,C 0O ,F : T) and Ci(R, C^, F, 7, r, T) 
such that for N > e ClT the following estimate is true 



Vte[o, T], Woo(M*),/2v(*))< 






(2.6) eq:thm3 



t m: cutoff 
result has also an interest for numerical simulation because one obvious way to 

fulfill the hypothesis on the infinite norm of /jy is is to put particles initially on a mesh (with 

a grid length of JV^ 1 / 2 ^ in M. 2d ). In that case, the result is even valid with 7 = 1. 



thm: deter 



thm: cutoff 



3 From deterministic results (Theorem E and 3) to 
propagation of chaos. 

Mun: deter 

The assumptions made in Theorem [B] may seem a little bit strange, but they are in some sense 
generic, when the initial positions and speed are choosen with the law (J )®^. Therefore, to 
prove Theorem fF ~ 



form Theorem 



we need to 



• Obtain a bound on the W\ distance between f{t) and iW(£), which are two solutions 
of the Vlasov equation. 

• Estimate the probability that empirical measure chosen with the law (f°)® N , do not 
satisfy the conditions i) and ii) of the deterministic theorem p , and are far away from f° 
in W\ distance (the last conditions is important for the previous point on the distance 
between / and f^). 

For these two points, we will use known results detailed in the next two sections. After that, 
a good choice of the parameter 7 and r will allow us to conclude the proof. 



3.1 Stability around solution of the Vlasov equation. 

i .opmOfi 
Loe06] for Qa = A — 1, but its proof may be adapted to our 



less singular case (The adaptation is done in |Hau09j in the Vortex case) 

Loeper Proposition 1 (From Loeper). If fi and fi are two solutions of Vlasov Poisson equations 
with different kernel K 1 and K 2 both satisfying a (S a )- condition, with a < d — 1, then 



d 
dt 



Wi(/i(t), / 2 (*)) < Cmaxdlpilloo, l|p 2 ||oo) Wx(fx{t), f 2 (t)) + C||pi|U||Xi - K, 



21 



The bound on the density may be obtained in our case with the argument of Pfaffelmoser 
for solution witfLcornpact support |Pfa92j . It is even simpler for a < 1 as it is explained in 
the appendix of I HJO'/jT - 

Using that theorem in the case without cut-off (K\ = K 2 = F) , with a < 1 (for d > 3) and 
a ll/jvll compactly supported, we obtain that there exists a constant C depending on F, an 
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uniform bounds on Jjie infinite norms of the /at and the size of their supports (denoted C^ 
and R in Theorem p ), such that 



W 1 (f(t)J N (t)) < e c ^W x {fj%) < e Cot (W.if^^ + N-^) 



(3.1) eq:Loeper 



3.2 Estimates in probability on the initial distribution. 

Deviations on the infinite norm of the smoothed ..empirical distribution fjj. The 

. prpp : larffsdev J 

precise result we need is given by the proposition [bjm the Appendix. It tells us that if the 
approximating kernel is < 



lr_i Ii2d, then 



P(||/£||oo>2 



l+2d 



,) < C 2 N^e 



-dN 1 —* 



with C 2 = (2R° + 2) d , R° the size of the support of /, and d = (2 In 2 - 1) 2 n 

We would like to mention that we were first aw^arenf the possibility of getting such estimates 
in a paper of Bolley, Guillin and Villani |BGV07[, where the authors obtain quantitative 
concentration inequality for \\f N — /||oo i n infinite norm under the additionnal assumption 
that f° and are Lipschitz. Unfortunately, they cannot be used in our setting because 
they would require a too large smoothing paramater. Gao obtain in l |Gao03] precise large 
deviations estimates for \\f N — /||oo, but as usual with large deviations estimates, they are 
only asymptotics and therefore less convenient for our problem, which is why we reworked 
them here. 



Deviations for the minimal inter-particle distance. It may be proved with simple 
arguments that the scale r] m is almost surely larger^ than N~ l l d when f° E L°°. A precise 
result is stated in the Proposition below, proved in |Hau09j : 

Proposition 2. There exists a constant c 2 d depending only on the dimension such that if 
f° E L°°(R 2d ), then 

We point out that this is not a large deviation result (the inequalities are in the wrong 
direction). It is that .condition that prevents us from obtaining a "large deviation" type 



that , , „ 

result in Theorem it (contrarily to the cut-off case of Theorem E ). In fact, the only bound it 
provides on the "bad" set is 



P \d N (Z) < 
With the notation of Theorem E 



I 



N l / d 

■ im: deter 



< 1 



e-^^ll- 1 * < c 2d \\fU d . 



it comes that if s = 7^ 



1 > then 



P (d N (Z) < e 1+r ) = P ( d N (Z) < 



jy-s/d 

N l / d 



< CMll/loo-nr*- 



(3.2) 



dN 
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Deviations for the W\ MKW distance. Peyre has obtained in IJPeypthe following result 

probaint Proposition 3 (Peyre). Assume that f° is a compactly supported measure on M. 2d . If d > 2, 
and the empirical measures (x° N are chosen according to the law {f Q )® N , then there exists an 
explicit constant Ld, depending only on the size of the support, such that 

P (wM,f) > ^)) < e*'"^ . (3.3) 

1 — 7 

It is of particular interest to us when L = N^~ , in which case it maybe rewritten 

^{W 1 {fi° N ,f°)>e)<Ce- N ^ L , with C = e Ld . (3.4) 

3.3 Conclusion 

■ tfchm:prob 

Now take the assumptions of Theorem IT It means that we assume that F satisfies a (S a ) 
condition for a < 1 and /° G L°° for d > 3. We chose 

f2 + 2a \ f 2 1 * d - l \ 

7 G I — ,1 , and rG --l,r = — — , 

\d + a J \7 l + aj 

the condition on 7 ensuring that the second interval is non empty. We also define 

1 + r ( d — 7 

*:=7— "1>0, A = min (^1 - 7, -^- 

Denote by u)±, u 2 the sets of initial conditions s.t. respectively (i), and (ii) (with the constant 
^ = 2 l+2d \\f \\ oo ) of Theorem gjhold and u 3 s.t. W^nJ ) < ^y. 

oji := {Z s.t. d N {Z) > s 1+r }, uj 2 := {Z s.t. ||/° l^ < 2 1+2d ||/ || oo } 
u; 3 :={^s.t. WM,f°)<e} 
By the results stated in the previous section, one knows that 

F(ul) < CN~ S , P(a£) < CN^e~ ClNl '\ P(a$ < Ce~ N ^ . 
Denote ui = iO\ fl uj 2 fl uj 3 . Hence \u c \ < \u>1\ + ju^l + 1^3 1 and for N large enough 

P(w c ) <CN' s + CN^e- m ' m(1 ' Cl)N ' X <CN- S (3.5) I boundomega 



■ fehm: deter 

If the initial conditions belong to ui then one may apply Theorem p and get on [0, T] 

ppot 
W X {f N ,H N ) < WooifN,^) < Nl/{2d) ■ 

J eq :Loeper 

Now apply the stability around solution of Vlasov equation given by (371) and get 



Wi(/, In) < W 1 (f°, f° N ) e°« < ^y e<*«. 

The factor 2 comes from the fact that W t (f° ', /&) < Wi(f°,n° N ) + Wi(jjP n , In)- We conclude 
that 

WiU^n) < j^- d e Cot , 

which proves that 

/ ^e Cot \ 

P(w) < P (W G [0, T], W 1 (f, f N ) < j^j-A . 

J houiidomega Ithm:prob 

The bound |B.5| then gives Theorem |[XJ 
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It hm: cutoff It hm : probcutof f 

3.4 From Theorem (Uto Theorem 

■ fchm : probcutof f J thm: cutoff 

In the cut-off case, one can derive Theorem [B]from Theorem g in the same manner. As we 
do not use the minimal distance in that case, the proof is simpler in the case a < d — 1 and 
we get a stronger convergence result. The only difference is that we shall use the Theorem^ 
with K\ — F and K 2 = -Fjv, so that an error term appears. But that error term is bounded 
by 

CWpflUlF - F N \U < Ce d - a < CWoo(t) 



eper 



for any t so that the proof is unchanged. In fact, with the same A, we obtain since the set 
lo\ is now useless that 



P 3te[o,T], W x (Mt),f(t))> 



4gCbi 



JV7/(2rf) J 
thm: deter thm: cutoff 



< CN~<e~ ClN * 



4 Proof of Theorem ||3j and 3 



4.1 Definition of the transport 

We try now to compare the the dynamics of jin and /at, two distributions which have a 
compact support. For that, we choose an optimal transport T^L=„TjL) from f% to n° N for 
the infinite MKW distance (See the remark after Definition ITF. The existence of such a 
transport is ensured by |CDPJ08bf. T° is defined on the support of f^, which is included 
in {\z\ < R } (the size of the support). Since /^ is an e-enlargement of [i° N , it is clear that 

Woo(/&,0<e. 

We also denote by Z? = (X*, V*) the smooth flow associated to /at and by Z N = (X N , V N ) 

the flow of the N particles system (with the convention Z(t, s) transport from time s to time 

t). A simple way to get a transport of /^(t) on /iAr(t) is to transport along the flows the map 

T°, i.e. to define 

T f = Z N (t,0)oT°oZ f (0,t), and T* = (T^T*) 

We use the following notation, for a test-" particle" of the continuous system at the position 
Zt = (xt,v t ) at time t, z s = (x s ,v s ) will be its position at time s for s G [t — T,t]. Precisely 



Z f (s,t, 



zt, 



Since f^ is the solution of a transport equation, we have /W(£, z t ) = /jv(s, z s ). And since the 
vector-field of that transport equation is divergence free 



/ $(z) f N (s, z)dz= f $(Z f (s, t, z)) f N (t, z) dz 



$(z s ) f N (t,z t )dz t . 



Finally let us remark that the f^ are solutions to the (continuous) Vlasov equations with an 
initial L°° norm and support that are uniformly bounded in N. Therefore this remains true 
uniformly in iV for any finite time. In particular there exists a constant C independent of iV 
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such that for any t G [0, T] 

IIMt,.,.)lloo<c, \\f N (t, 



|£?U(t) := ||S(t,.)||oo<sup / |F(x-y)|/jv(t,y,«;)dj/dt«<C 

X J 

\VE(t,x)\< [ \VF(x -y)\ f N (t,y,w)dydw\ < C 



(4.1) 



boundfN 



suwf N (t,.,.)eB(0,R(t)), R(t)<C, 



as of course R(t) < R° + J \\E(s, .)||oo ds. This is always true for a < 1. In dimension d < 3 
it remains true for a < d — l^and even a = d — 1. In fact, all that estimates where central in 
the work of see Pfaffelmoser l |Pfa92] about existence, and uniqueness of compactly supported 
solution of Vlasov-Poisson equation (See also l [Hor93] for a result with improved bounds). 
The proofs can be adapted to our simpler cases (See the Appendix of |HJ0'7jror the case 
a < 1). 

In dimension d > 3 and for attractive forces with 1 < a < d — 1, there can be a blow-up in 
finite time (for a larger than a critical value depending on the dimension). In that case, we 
simply restrict ourselves to a time interval on which this does not occur. 

In what follows, the final time T is fixed and independent of N. For simplicity, C will denote 
a generic universal constant, which may actually depend on T, the size of the initial support, 
the infinite norms of the /at... But those constants are always independent of iV as in (S~T) 



4.2 The quantities to control 

We will not be able to control the infinite norm of the field (and its derivative) created by 
the empirical distribution, but only a small temporal average of this norm. For this, we 
introduce in the case without cut-off a small time step r = e r for some r' > r and close to r 
(the precise condition will appear later). 
In the case with cut-off where r and r' are useless, the time step will by r = e. 

Before going on, we define some important notations. 

• The MKW infinite distance between //jv(£) an d /(£)• 

We of course wish to bound Wodt) := sup 0<s<t W 00 (//jv'(s), /jv(s)) (note that Wca is 
hence automatically non decreasing). For the transport introduced before, one has 

Woo(t) < sup sup \T s (x s ,v s ) - (x s ,v s )\. 

s<t (x s ,v 3 )£supp f N (s,.,.) 

In fact, we will provide bound for the quantity of the right hand side. Our result maybe 
stated for that quantity, rather than the infinite MKW distance. It is a little stronger, 
since it means that rearrangement in the transport are not necessary to keep the MKW 
distance bounded. The transport chosen at time t = is preserved during the time. 

• The support of fj,^ 

We shall also need a uniform control on the support in position and velocity of the 
empirical distributions : 

R N (t) = m a x\(X t (t), Vi(t))\. 
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The infinite norm l-E^loo of the time averaged discrete force field. 

We also define the average of the discrete force field over small time intervals of length 
r (the dependence on t is implicit) 

1 I* 



| .Eat |oo = sup- / \E N (Xi(s))\ds 



The infinite norm iV^-Eloo of the time averaged discrete derivative of the 
force field. 

We also define a version of the infinite norm of its averaged derivative 
7iVci| 1 A* \E N (X t (s))-E N (X 3 (s))\ds 



|V £|oo = sup - / — . . ds. 

<# r J t _ T \Xi(s) - Xj{s)\ + e^ l +' > 

For both E^ and V E, we use the convention that when the interval of integration 
contain (for t < r), the integrand is null on the left side. The control on that term is 
useless in the cut-off case. 

The minimal distance in phase space d^ 

■l eg :d min J sec :main 

which has already be defined by the equation (12. if in the Section 12) 



Two useful integrals I a (t,z t ) and J a+1 (t,z t ) 
Finally the technical computations involve 



I a (t,z t ,z t ) = - [ \F(T s x {z 8 ) -T*{z a )) - F(x s - x.)\ds, 

T Jt-T 

which controls the difference of the two force fields at two point related by the "optimal" 
transport. Defining a second kernel as 



K F = min 



L L" 



1 



J.I1 + Q ' pl+5"' \rn\ 



we define a second useful quantity 



i-t 
J a+1 (t,zt,zt) =-- - / K e {\T s x {z s )-T s x {z s )\ds 

ft 

K^X^-X^ds, 



t-r 



if i and j is the indices such that Zi(t) = T*(^) and Zj(t) = T'(zi). J a +i will be useful 
to control the discrete derivative of the field, and is thus useless in the cut-off case. 

All previous quantities are relatively easily bounded by I a and J a +i- Those last two will not 
be bounded by direct calculation on the discrete system, but we will compare them to similar 
ones for the continuous system, paying for that in terms of the distance between //jv(£) and 
/(£). That strategy is interesting because the integrals are easier to manipulate than the 
discrete sums. 
We summarize the first easy bounds in the following 



Mm: deter 

Proposition 4. Under the assumptions of Theorem^ one has for some constant C uniform 

in N 

(i) R N (t) < Woo(t) + R(t) < W^(t) + C, 



{it) Wooit) < W^t - t) + C t swp / I a (t,Zt,Z t )dZt, 

zt J\zt\<R(t) 



propeasy 



(Hi) |V E\oo < C sup / J a +i(t, z t , z t ) dz t . 

zt J\z t \<R{t) 

(iv) d N (t) + e 1+r ' > [d N (t - r) + £ i+r'] e -(i+|v^Mt)). 



Note that the control on Rx(t) is simple enough that it will actually be used implicitly in 
the rest many times, and that the iv) is a simple consequence of the Hi). In fact, in that 
proposition the crucial estimates are the ii) and Hi). 

Remark also that in the case of very singular interaction force (a > 1) with cut-off - in short 
(S^) conditions - the control on minimal distance d^ and therefore the control on iV^-Eloo 
are useless, so that the only interesting inequality is the second one. 

i propeasv 



ror 

4.3 Proof of Prop, gf 

Let us start with (i). Simply write 



R N (t) = sup \T\z t )\ < sup \T t (z t )-z t \ + sup \z t \, 

z*esupp/jv(i,-) z t esupp/jv(i,-) z t esupp/jv(i,-) 



■ Ibou j . 

So indeed by (ffTT[ 



atmdfN 



R"(t) < W^it) + R(t) < W^it) + C. 



As for (ii), simply differentiate in time W^ to find 



W^-W^t-r) 



< sup 

zt 



[F(T:(z s )-T:(z s ))-F(x s -x s )]ds 



t-T 



f N (z t )dz t . 



Since /at is uniformly bounded in L°° and compactly supported in B(0,R(t)), one gets by 
the definition of I a 



vM*)-w°o(*-t) 



< ||/°||oo sup / I a (t,z t ,z t )dz t , 

zt J\z t \<R(t) 



which is exactly (ii). 

Concerning iV^-Ejoo in (Hi), noting that 



' \E N {Xi{s)) - E N {Xj{s))\ 1_ 
tT [X^s) - Xj(s)\ + e^' N 



E 



\F(X l ( S )-X k (s))~F(X,(s)-X k ( S ))\ 



+ 



1 f* \F(X t (s) - X 3 (s)) - F(X ] (s) - X^s)^ 



N 



t-T 



IXAs) - Xi(s)\ + e^' 



ds. 



ds 
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i aa : . Ualpna 
By the assumption (TTo), one has that 



So 



l'^--^l^l^p + ^H J ' "I 



\F(X t (s) - X k (s)) - F(Xj(s) - X k (s))\ < C C 



\X i { 8 )-X i { 8 )\+£^ ~ \Xi{s) -X k (s)\i+" \X,(s) -X fe ( s )|i+°' 

and that bound is also true for the remaining term where _k = i or ?', if we delete the undefined 

Jeq-^C alpha J ' 

term in the sum. One also obviously has, still by (IXToJ 



iFiX^s) - X k (s)) - F(X 3 (s) - X k (s))\ < _ C 



\X t (s)-X 3 (s)\+e^ r ' -e^lXiW-X^s) 

C 

+ 



ei+r'\ Xj (s)-X k (s) 
Therefore by the definition of K e 

iF^s) - X k (s)) - FiX^s) - X k (s))\ 



\X t (s)-X,(s)\+e^' 
Summing up, this implies that 



< K^X^s) - X k (s)) + K^Xjis) - X k (s)). 



jN E\ 

J -'\oo 

*# VT„,_, 

kjtt 



T J t-r N U 

+ \f ^J2 K e(X 3 (s)-X k (s))ds). 



Transforming the sum into integral thank to the transport, we get exactly the bound (Hi) 
involving J a +i- 

Finally for d^(t), consider any i ^ j, then obviously 

fj(Ms) - Xj(s), Vi{s) - Vj(s))\ > -mis) - V,(s)\ - \E N {Xi{s)) - ^(X,-(s))|. 
Simply write 

\e h (M»)) - WWH < ^^r^f^P (l*M - *,(.)| + ^ 1+ "> 

to obtain that 

d UY(\ Y(\V(\ T/^M*s fi \ E N(X 1 (S)) - E N (X,( S ))\\ 

(KX^-X^V^-VjWl+e 1 ^'). 
Integrating this inequality and taking the minimum, we get 

A U\i l+r'^M U ^ 1+r'^-f ( C ^^(s)) - E N (Xj(s))\ A 

d N (t)+e + >(d N (t-r)+e + )miexp(-T- — — J ds 

> [d N ( t _ r ) + £ l+r'] exp -r(l+|V^ E U(t)) 
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lem:boundI 



boundlK 



boundlcut 



4.4 The bound for I a and J a +\ 

S opeasy 
it remains to bound the two 

integrals involving I a and J a . It is done with the following lemmas 

Lemma 1. Assume that F satisfies an (S a )- condition with a < 1, and that r is small enough 
such that 

Ct (1+ \V N E\^{t)) (Woo(t) + r) < d;v(t). (4.2) 

Then one has the following bounds, uniform in z t 

I a (t, zt, z t ) dz t < C [W^t) + (Woo(t) + r) d r- a + (W^t) + r) 2d (d N (t))- a r- a ] . 
\z t \<R(t) 



J a+1 (t, z t , z t ) dz t < C (1 + (W^t) + r) d e^ l+ ^ r~ a 

\z t <R(t) 

+ (Woo(*) + r) 2d e^ 1+ ^ r~ a (d N (t))- a ). 



In the cut-off case where the interaction force satisfy a (S^) condition (we recall that it means 
that the cut-off is of size N~ m = e m with rh = —m), we only need to bound the integral of 
I a , with the result 

Lemma 2. Assume that 1 < a < d — 1, and that F satisfies a (S^) condition, one as the 
following bound, uniform in z t 

[ I a (t, z t , zt) dzt < C (W^t) + (W^t) + r )<V-i £ -(i-") + (WUt) + T) 2d e- ma ) . (4.3) 

J\z t \<R{t) 

with the convention (if a = 1) that (e m )° = \ ln(e m )| FJ 

The proofs with or without cut-off follow the same line and we will prove the above lemmas 
at the same time. We begin by an explanation of the sketch of the proof, and then perform 
the technical calculation. 

4.4.1 Rough sketch of the proof 

The point z t = (x t , v t ) is considered fixed through all this subsection (as the integration 
is carried over z t = (x t ,v t )). Accordingly we decompose the integration in z t over several 
domains. First 

A t = {z t | \x t - x t \ > 4Woo(i) + 2r(\v t - v t \ + r|£|oo(*)) }■ 

This set consist of points z t such that x s and T£(z s ) are sufficiently far away from x s on the 
whole interval [t — r,t], so that they will not see the singularity of the force. The bound over 
this domain will be obtained using traditional estimates for convolutions. 

x That convention may be justified by the fact that it implies a very simple algebra (x 1-Q )' ~ x~ a even if 
a = 1. It allows us to give an unique formula rather than three different cases. 
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cond:leml 



boundlKcut 



One part of the integral can be estimated easily on A c t (the part corresponding to the flow 
of the regular solution f^ to the Vlasov equation). For the other part it is necessary to 
decompose further. The next domain is 

B t = A c t f){z t | \v t - v t \ > 4Woo(t) + 4r|£| 00 |(t)}. 

This contains all particles z t that are close to z t in position (i.e. x t close to x t ), but with 
enough relative velocity not to interact too much. The small average in time will be useful 
in that part, as the two particles remains close only a small amount of time. 
The last part is of course the remainder 

C t = {A t U B t ) c . 

This is a small set, but where the particles remains close together a relatively long time. 
Here, we are forced to deal with the corresponding term at the discrete level of the particles. 
This is the only term which requires the minimal distance in phase space; and the only term 
for which we need a time step r small enough as per the assumption in Lemma pi 
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Figure 1: The partition of the phase space. 



4.4.2 Step 1: Estimate over A t 

If ZfgAf, we have for s G [t — r, t] 

\x s -x s \ > \x t - x t \ - (t - s)\v t - v t \ - it - s) 2 |£|oo(£) > 
\T a x {z s )-T s x {z s )\ > \x s -x s \-2W OD (s)> \ Xt ~ Xt V 



\%t ~ x t\ 



(4.4) 
(4.5) 



For I a , we use the direct bound for z t G A t 

\f(t:(z s )-t:( Zs ))-f(x s -x s )\ < 



c 



x t - x t 



1+a 



(\T*(z s )-x s \ + \T*(z s )-x s \) 



< 



Woo(s) < C |1+tt Woo(t), 

\xt — xA + 



\x t — x t \ 1+a 
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cylinder 



and obtain by integration on [t — r, t] 

\ x t — Xt\ + 

Then integrating in Zt we may get since a + 1 < 2 < d 

dz t 



J At J A 



A t \x t -x t \ 1+a (4.6) [kt_ 

< CR{t) 2d - 1 - a W 00 {t) < CW^t). 



FAR 



For J a +i, we have used (KoTon the set A t the bound 

\K £ (T°(z t )-Ts(z t ))\<- ( 






Integrating with respect to time and z t we get since 1 + a < d. 

Ja+i{t, %&, z t ) dz t < C 



A t "'* J At \x t -x t \ 1+a (4.7) [kt2_ 

< CR(t) d - l - a <C. 

For the cut-off case, the estimation on I a for this step is unchanged. 

4.4.3 Step 1' : Estimate over A\ for the continuous part of I a and J a +i- 

For the remaining term in /„, we use the rude bound 

\F(T:(z s ) - T°{z s )) - F(x s - x.)\ < \F(T:(z s ) - T*(z s ))\ + \F(x s - x s )\. 

The term involving T s is complicated and requires the additional decompositions. It will be 
treated in the next sections. The other term is simply bounded by 

1 \^,- m , , 1 /"* f Cdz t , 
\F(x s — x s )\ds dzt < — I / -p r^ ds 

z t ^A c t ' Jt-r T Jt-r Jzt<=A c t \ x s ~ x s\ a 

< l -f I F^ *- 

T Jt-r Jz B eZf(s,t,Ac) \ x s %s\ 

From the bound R(t) and l-E'loo(^) we see that 

\A c t \ < CRitfiW^it) + r) d < CiW^it) + r) d , 

where | • | denote the Lebesgue measure. Since the flow Zf is measure preserving, the measure 
of the set Zf(s,t,Af) satisfies the same bound. This set is also included in [— R(t), R(t)] 2d 
(if R is increasing, a property that we may assume). We use the above lemma which implies 
that above all the set Z(s, t, Af), the integral reaches is maximum when the set is a cylinder 

Lemma 3. Let Q G B(0, R) G W 1 . Let P be a projection from M. n to M. m with m < n. Then 
for any a < m 

I U,dj j~>a(n—m)/m | r~) 1 1 — aim 

Jn l ri l 
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icylinder 
Proof of Lemma pL We can freely assume that Px = (x±, . . . ,x m ). Now maximize the 

integral 

-i tAj \AJtAJ 

over all sets u C M n satisfying u> C B(0,R) and \u>\ = \Q\. It is clear that the maximum is 
obtained by concentrating as much as possible u near Px = 0, i.e. with a cylinder of the 
form B m (0, r) x i? n _ m (0, R) where B^ denotes the k-dimensional ball. Since \u\ = \Q\ we have 
r = \il\ 1 / m R 1 ~ n ^ m . The integral can now be computed explicitly and gives the lemma. □ 

Applying the lemma, we get 

- / \F(x a - x s )\ dz t ds < CiW^t) + r) d - a . (4.8) 

JLem; boundl 

That term do not appear in Lemma [LL] since it is strictly smaller than the bound of the 

remaining term (involving T), as we shall see in the next section. 

The same bound is valid for J a +i since a + 1 < d. 

For the cut-off case, the same bound is valid for I a since a < d — 1 < d (The cut-off cannot 

in fact help to provide a better bound for this term), and we do not need the estimate on 

Ja+l- 

At this point, the remaining term to bound in I a is only 

L ? L m ° ] - r - fe))i ds s c L l r L pctf ds - (4 - 9) 

For J a+ \ one may bound the integral of the continuous part on A c t in a similar manner. As 
K e < i+X ia ; the remainder can be controlled by (IO) 



J a+1 (t,z t ,z t )dz t < C(W 00 (t)+r)- a 4 — T I - I ]Ts( -r Ts(z]]a d^ ( I 10) 

A1 £ Jz t eA c t - T Jt-T \ 1 x\ Z s) ~ 1 x\ Z s)\ 



vd-« , c f i r ^ t 



,\ ea _;jti 

Therefore in the next sections we focus on giving a bound for (47y JT 



4.4.4 Step 2: Estimate over B t 

We recall the definition of B t 

\x t - x t \ < 4W 0O (t) + 2r(\v t - v t \ + rl^looCt)) 



Bt \ Zt S,t - |^-^|>4W oo (0+4r| J B| oo |(t) 
If Zj G 5t, we have, as for A t , for s G [t — r, t] 

\v 9 -v,-v t + v t \ < 2r\E 00 \(t)< \ Vt ~ Vt \ (4.11) 

\T s v {z s )-T°{z s )-v t + v t \ < \v s -v s -v t + v t \ + 2W 0O (s)\ < ^ ~ ^ . (4.12) 

This means that the particles involved are close to each others (in the positions variables), 
but with a sufficiently large relative velocity, so that they do not interact a lot on the interval 

[t-T,t]. 
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Ic 



eq: remain 



remainderJ 



m 



First we introduce a notation for the term of ( 



I bc (t, z t , z t ) dz t , with I bc (t, z t , z t ) = I bc (t, i,j) ■= - 

zt£B t T 



where (i,j) are s.t. T£(z s ) = X^s), T°(z s ) = Xj(s). 
For Zt G B t , define for s G [t — r, t] 



t- r \T£{z a )-T»(z,)\« 



ds, 



(4.13) eq:Ibc 



(«) = CS&) - t:(z s )) 



Vt -vt 

\v t - v t \ 



Note that \<f>{s)\ < \T*{z 9 ) - T°(z s )\ and that 



<f>'(s) = (T*(z s ) - T*(z s )) 



vt-vt 
\v t - v t \ 



fa - v t \ + (T'(z.) - T„ s (* s ) - ( Vt - Uj)) • : > 



where we have used (14.11 |T" Therefore if attains its minimum on [t — r, t] at sq, then 



|T s (^ s )-T s (^)|>|t- So | 



ft -fd 



(4.14) eq:dispB 



Using this directly gives, as a < 1 



C 
\I bc (t,z t ,z t )\ < —\v t -v t \ 

T 



ds 



t- T \s — s | 



Now integrating 



\Ibc(t, zt, zt)\ dzt <Ct 

z t eB t J A c t \ v t ~ v t\ 



< C T- a \v t - V t \- a . 



dz t 



(4.15) 



vtrick 



<Cr- a (W 00 (t) + r) d (R(t)) d - a 1 
by using the fact that B t C B(0, C[Woo(t) +r]) x B(0, R(t)). In conclusion 

|/ 6c (i, z t , zt)\ dz t < C r- a (Woo(t) + r) d . 



(4.16) 

zteBt 

Jvtripk 
With the cut-off where a > 1, the reasoning follows the same line up to the bound (14.15 P 

J vtripg l i r 

which relies on the hypothesis a < 1. ( 14715" pis replaced by 



Bt 



C /"* 

|/te(t,^,zt)| < — / 

7" it-r 



<is 



(\s - s \\v t -v t \ +4e T 



< 



c n 



T 



ds 



(s + Ae" 



< 



Ce fh{l-a) 



When a — 1, the previous calculation gives 

C f T 

\I bc (t,z t ,z t )\ < — In 1 + — 

r V 2\v t -v t \e 



< 



Cln(e-™) C(e" 



m 

rh\0 



)<7ln(l+^) 



r 
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In the first line, we used the bound t^ — , < ,„, 1 u ^ which is implied by the definition 

' \Vt— Vt\ — |-E|oo(*)) l J 

of B t and then we assumed that \E\ 00 (t) > 1. In the second line, we use the convention 

(e*)°= |ln(e*)|. 

In both cases, integrating that bound (which do not depend on Vt) over B t , we get the 

estimate 

/ \I bc (t,z t ,z t )\dz t <C(W 00 (i)+T) d T- 1 e*( 1 - a >. (4.17) 

Jz t eB t 

4.4.5 Step 3: Estimate over C t 

We recall the definition of C t 

r _ j , \x t -x t \< 4W 0O (t) + 2r(\v t - v t \ + r|£|oo(*)) 
1 \ Zt |v t -v t |<4W 0O (t)+4r|SU(t) 

First remark that C t C {\z t \ < C(W 00 (t) + r)}, so that its volume is bounded by CiW^it) + 
r) 2d . From the previous steps, it only remains to bound 

hc(t, Zt, z t ) dz t . 
z t ec t 

We begin by the cut-off case, which is the simpler one. In that case, one simply bounded 
he < C e~ ma which implies 

I bc (t, z t , z t ) dz t < CiW^t) + r) 2d e- ma . (4.18) 

Zt€Ct 

It remains the case without cut-off. We denote C t = {j \ 3z t G C t , s.t. Zj(t) = T*(^)}, and 
transform the integral on Ct in a discrete sum 

\ I bc (t,z t ,z t )dz t = y2a i:j I Nc (t,i,j) with I Nc (t,i,j) = - —— - — * ds, 

JztGCt . ?. T Jt-nr l A iV s J — A jV s JI 

3&C t 

where i is the number of the particle associated to z t (^(zt) = Zi(t)) and 

ciij = \{z t e C t , T l (z t ) = Zj(t)}\, so that ^ a„ = \C t \. 

jec t 

To bound In c over Ct, we do another decomposition in j. Define 

JX t = ^eCt,\x,{t)-xm> dj Y^y 

d N (t) 



Bt cut off 



Ct cut off 



JV t = <jeC t , \Xj(t) - Xi(t)\ < \Vj(t) - Vi(t)\ and \V 3 (t) - V t (t)\ > 

By the definition of the minimal distance in phase space d^it), one has that Ct = JX t U JVt- 
Since 

^(z^-Ztl^Wooit), 

one has by the definition of C t and C t that for all j e C t , \Zj(t) — Zi(t)\ < C {W^it) + r). 
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Let us start with the bound over JX t . If j G JX t , one has that 

|X,-( S ) -^(s)! > |X,(t) -X,(t)| - f \V,{u)-V % {u)\du. 

J s 

On the other hand, for u G [s,t], 

\Vj{u) - Vi(u)\ < 2W 0O (t) + \v u - v u \ < 2(W o0 (t) + rl^U) + \v t - v t \ < C{W^{t) + r). 
Therefore assuming that with that constant C 

CriW^ + ^KdNit)/^ (4.19) 

we have that for any sG [t — r, t], \Xj(s) — X{(s)\ > djy(t)/4. Consequently for any j G JX t 



W*,i,J)<C[M*)]" a - 



(4.20) 



assJX 



boundJX 



For j G JV^, we write 



1(^(8)-^))- (^(t) -Vi(t))l< / IM*i(«))--M*i(«))l*'- 



Note that 



|X i (8)-X i (s)|<|X i (t)-X i (t)|+ / |^(w)-^(«)|d« 



< C(Woo(t) + r) + 2 / (WooH + #(«)) rfw ( U ' ' 

<C(^oo(*)+t). 

Hence we get for s G [t — r, t] 

I \E N (X s (u)) - ^(^(«))| du < Ct\V n E\ oc (W^t) + r + e 1+r '). 

J s 

Note that the constant C still does not depend on r = e l+r . Therefore provided that with 
the previous constant C 



xclose 



Ct IV^U (W^t) + r) < d N (t)/A, 



(4.22) 



assJV 



one has that 



\Vjis) - V t {s) - (V#) - Vi{t))\ < d N (t)/A and also \V t (s) - Vj{s)\ > 



d N (t) 



J gq :dispB 

As in the step for B t (See equation (4.14)) this implies the dispersion estimate 



\Xj(s) — Xi(s)\ > \s — so| djsr(t)/4: for some sq G [t — r, t]. As a consequence for j G JV t , 



C 

I Nc (t,i,j) < —(d N (t))~ 

T 

. --UJX . Ibqund JV 

Summing (4.20p and (|4.23rrone gets 



ds 



-T \ S ~ S 0| 



<Cr-° (<**(*))" 



(4.23) 



boundJV 



J2 aijlNcit, i,j) <C\C t \ ((d N (t)y a + r- a (d N (t))- a ) . 
jec t 

Coming back to he, using the bound on the volume of \Ct\ and keeping only the largest term 



of the sum 



I bc (t, zu zt) dz t < C (W^t) + T) M r~ a (d N (t)) 



(4.24) 



Ct 



c t 
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4.4.6 Conclusion of the proof of Lemmas 






tofiffiKdlcut 



Assumptions ( 4 l7 IJ^ r and (14.22 r are ensu r e^ by the hypot hesis of the lemma. Summing up 
(lO) for I a or (577 1 for J a+ \, with (478), (4716) and (4724), we indeed find the conclusion of 



the first lemma. 

J 



of 
EtetLtoff 



ILL-. XLc-. JBtcujt 

Inthe S" case, no assumption is needed, and summing up the bounds ( |4T6[ ), ( |478[ ), ( [4717 ) 
(l4.18|), we obtain the second lemma. 



Jthm: deter 

4.5 Conclusion of the proof of Theorem |S] (without cut-off) 



In this subsection, in order to make the argument clearer, we number explicitly; the constants. 
Let us summarize the important information of Prop. |4| and Lemma |i| Let us also rescale 
the interested quantities s.t. all may be of order 1 

eW^t) = Woo(t), e 1+r d N (t) = d N (t). 

We also assume that W^ and iV^-Eloo are non decreasing functions of t, and that d^ is a 
non increasing function of t. In fact the bound proved before are also valid for sup s<4 W^,, 

sup s<t iV^-Bjoo, and inf s<t d^^s). With that convention Woo(t) > WooiO) > 1. By assump- 
— i fehmTdeter ~ 

tion m in Theorem B also note that d N (0) > 1. _ 

w ,l—i iv \/— [ Hrem : boundl 

Recalling r = e r (with r > r > 1), the condition of Lemma 111 after rescaling reads 



d e r r (1 + |V JV £U(t)) Wooit) < d N (t), V t E [0, to]. 



JLem : boundl 

In Lemma |ll| we proveol that there exists some constants Cq and C\ independent of A" (and 



hence e), such that if (4.25) is satisfied, then for any t e [0, to 

Woo(t) < W^t -r) + C e r ' (w^t) + e Al Wtit) + £: A2 W 2 J{t) ,/-"(/)) . 

\V N E\ OQ {t) <C 2 (l + e x * Wtit) + e x * W 2 J(t) d N a (t)) 



T 



£ r'-n e -r(l+\^ N E\ 00 (t)) 



d N (t) + e r - T >[d N (t 
where e appear four times with four different exponents A,,i 



, 4 defined by 



Ai = d — 1 — ar', 

As = d — 1 — r' — a r' 



A, 



A 2 = 2d -I- a(l + r' + r), 

-Id — \ — r — a(l + r' + r). 



To propagate uniform bounds as e — > and A" — >■ oo, we need all A« to be positive. As 
r, r' > 0, it is clear that Ai > A3 and A2 > A4. Thus we need only check A3 > and A4 > 0. 
As r' > r, it is sufficient to have 



r' < 



1 + a' 



and r' < 



2d 



a 



l + 2a 



Note that a straightforward calculation shows that 



d 



2d 



a 



a 



d 



1 + a l + 2a (l + a)(l + 2a) 



<0, 



(4.25) assumption 
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so that the first inequality is the stronger one. Thanks to the condition given in Theorem I 



im: deter 



d- X 

1+ 



-, so that if we choose any r' G (r, r*), the corresponding Aj are all positive. We 



fix a r' as above and denote A = minj(Aj). Then by a rough estimate 

WU*) < ^oo(t - t) + C r (^(t) + 2e x W 2 J(t) d~ N a (t) 
\V N EUt) < C 2 (l + 2e x W 2 J{t)d N a {t)) , 

d N (t) > [d N (0)+e r '- r }e-^ 1+ ^ NE ^^-e r '- r . 

■ [assumption 

If one has (14.25) and 



2e x W™(t)dJ?(t)<l, 
then using W^ > 1, we get W^it) < W^it - r) + 2C rW / " oo (t) so that 

W oo (t)<W oo (t-r)(l-2C r)-\ 

\V N EUt)<2C 2 , 

d N (t)> e -(^C 2 )t_ £ r>-r^ 

The last inequality implies d N (t) > I e " (1+2C2)t if 2e T -' r ' e^ 2C ' l)T < 1. That condition is 

fulfilled for e small enouetu i.e. N large enough (lnA^ < T). 

The first inequality in (Eggf; iterated gives W^t) < W / " O o(0)(l - 2C r)-r. If C r < |, then 



(4.28) 



we can use — ln(l — x) > 2x for x G [0, ^], and get 



Woo{t) < ^oc(0)e 



4C t 



To summarize 



^oo(t) < e 4Cof , 



V"E\ 00 (t)<2C 2 



(4.29) 



djv(*) > 2 e 



-(l+2C 2 )t 



At the discrete level of the particles, the dynamics is continuous in time, at least for initial 
conditions not leading to collisions., That set is of full measure for a < 1, and d > 2 (See 

H annzf ° assumption as_s. i J.'nption2 , v 

l |Hau04j ) . So as long as (14.25) and (14.27) are satisfied at t = 0, there exists a maximal time 
t G]0, T] (possibly t = T) such that they are satisfied on [0, to]. 

We show that for AMarge enough, i.e. e small enough, then one necessarily has to — T. 
i (114.291) on [0, T] which is the desired result. This is simple enough. 



Then we will have (4.29 ron [0, T] 



igh. By 



contradiction if to < T then 

C x e {r - r,) (1 + | V^UCto)) ^oo(to) = d N (to), or 4 e A W™(t ) d~ N a (t ) 

■Ire suit 



Until t , ( 



lolds. Therefore 

e A iy^(t )V(t ) <^2 Q e (a + (4d+2a)max{Co ' C2))40 < 1, 



■ [assumption 

for e small enough with respect to T and the Cj. This is the same for (14.25 F, 



C x e^-^ (l + l^EUt^W^to^ito) < 26^0,(1 + 2C 2 )e( 1+6 ^ Co > c ^ to < 1. 

■ fchm: deter 

Hence we obtain a contradiction and prove Theorem |B 



(4.26) roughestim 



(4.27) assumption 



presquegro 



result 
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Jthm: cutoff 

g (cut-o 



4.6 Conclusion of the proof of Theorem ||4j (cut-off case) 

In the cut-off case, using Lemma Eftogether with the inequality ii) of the Proposition Er 
may obtain 



opeasy 

we 



Woo(t) <W 00 (t-r) + CoW^it) [l + (W^t) + rf-'r-'e^ 1 -^ + (W^t) + r) 2 ^^"*"] . 

We again rescale the quantity W^it) = eW^t) and replace W^it) by W^it) + 1. Choosing 
in that case r = e, it comes for 1 < a < d — 1, 



WU*) < Woo(* - t) + CoWooCt) [l + e *-2-ft(a-D ly^-i(t) + £ 2*-i 



"™^- 1 (i) 



As in the previous section, we will get a good bound provided that the power of e appearing 
in parenthesis are positive. The two conditions read 



m < m* 



mm 



d - 2 2d - 1 



o 



o 



In that case, for iV large enough (with respect to e ), we get a control of the type 

jWoo(t) < 4C Woo(t), 

(but discrete in time) which gives the desired result. 

Remark 2. In t/ie cut-off case (and also in the case without cut-off'), it seems important to 
be able to say that the initial configurations Z we choose have a total energy close from the 
one of f° . Because, if the empirical distribution [i z N is close form f°, but has a different total 
energy, we would not expect that they do not remains close a very long time. 



: cutoff 

Fortunately, such a result is true and under the assumptions of Theorem [G] and gj the total 
energy of the empirical distributions are close from the total energy of f°. 
Unfortunately, we do not have a simple proof of this fact. But, it can be done using the 
argument for the proof of the deterministic theorems. First, the difference between the kinetic 
energy is easily controlled because our solutions are compactly supported and that there is no 
singularity there. Next, performing calculations very similar to the ones done in the proofs, 
we can control the difference between a small average in time of the potential energies, on 
the small interval of time [0, r] . Then, we control the average of the total energy, which is 
constant. 



p:largedev 



A Appendix : Large deviation on the infinite norm of 
In- 

Proposition 5. Assume that p is a probability on W 1 with support included in [— R°, R°] n and 
and bounded density f(x)dx. Let be a bounded cut-off function, with support in [— f, f] n 
and total mass one, and define the usual <p £ := ^0(-). For any configuration Z^ = (Zj)j<n 
we define 
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If e(N) = iV™ and the Zn are distributed according to f® N , then we have the explicit "large 
deviations" bound with c = (2L) n ||0|| oo and c = (2R° + 2) d L~ n . 



F(||/-|U > ^ll/IU) < c 7V7 e -(^n^ + i)(2L)»||/|Uiv 1 -,_ 
In particular, for <fi = l[-i/2,i/2] n and (5 = 2, we get 

P (Halloo > 2 1+n ||/|| 00 ) < (2W + 2) d We-^- l ™\~ Nl -\ 



(A.l) eq:largede 



(A. 2) eq:largede 



Proof. For any Z G 



l nN and z G R", we have 

jV 



/£(*) 



a7^> 



'-2 - Z; 






Ae 



11/. 



JVlloo 



< 



< 



Ne 



i=l ~ i=l 

#{is.t. |*-Z»|oo<¥} 



2 - Z,; 



v - SUp #{? S.t. |z - Z^ < -f }, 



where # stands for the cardinal (of a finite set). It remains to bound the supremum on all 

the cardinals. The first step will be to replace the sup on all the z G M. n by a supremum on 

a finite number of point. For this, we cover [— R°, R°] n by M cubes Ck of size Le, centered 

at the points (ck)k<M- The number M of square needed depends on A via e, and is bounded 

by 

"2(i?° + l) 



M < 



Le 



Next, for any z G M. , there exists a k < M such that \z — Ck\ < -§-■ This implies that 



sup #{i s.t. \z - Zilvc < -f} < sup #{i s.t. \c k - Z^ < Le} 



.rfcL 



fc<Af 



Now we denote by if^ := #{i s.t. |c£ — Zj|oo < Le}. H™ follows a binomial law B(N,p k ) 
with p,fc = f 2C f(z) dz, where 2Ck denotes the square with center c k , but size 2Le. Remark 
that 

Pk<p:=(2Le) n \\f\\ 00 . 

For any A, the exponential moments of H^ are therefore given and bounded by 



E(e XH ") = [l + (e x -l)p k ] 



N 



< [l + (e A -l)(2L£) 

< e< 



i jv 



j(e x -l)N(2Le) r 



Now for the supremum of the H^ 



E(e XsuPkH ") < E(e XH ?) + --- + (e XH M) 

< Me (eA " 1)Af(2Le)n|l/l|o ° 
'2(R° + V' n 



< 



Le 



(e A -l)/V(2L e )"||/|| c 



31 



Using finally Chebyshev's inequality, we get for any (3 > 



P(||^|U>/3(2Lr 



co U Wool 



< F(supH?>f3\\f\\ocN(2Le) n 

< E / e Asu Pfe ^fx -A/3||/|| 00 7V(2L £ )" 



< 



2(R° + 1) 
Le 



e A -l 



-\p)N(2Le) n \\f\\ c 



For > 1, the optimal A is In /3 and we get with c^ = (2L) 

"2(i?° + 1 



P(||/£||oo>/^ 



< 



Ls 



-(/31n/3-/3+l)Af(2Le)"||/|| 00 



With the scaling e(iV) = N », we get 

Remark finally that the choice of scale e(iV) = (In N)N~™ is also sufficient to get a probability 
vanishing faster than any inverse power. □ 
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